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Abstract 

Bilocal correlator in gluodynamics is evaluated at finite temperature in the framework 
of instanton gas model. It is demonstrated that vacuum correlation length decreases with a 
growing temperature. Obtained results are compared with lattice data for the bilocal correlator 
at finite temperature. Density of instantons and possible structures of nonperturbative vacuum 
are discussed. 

1 Introduction 

Nonperturbative fields play important role in the vacuum of QCD and gluodynamics, and determine 
many properties of the theory, in particular, they are responsible for the confinement, spontaneous 
chiral symmetry breaking (SCSB), and as a result, for the formation of observed hadron spectrum. 
A well developed and successful in the description of many properties of QCD model is the instanton 
liquid model, suggested by Shuryak T and Diakonov and Petrov The main nonperturbative 
fields in this picture are well separated instantons and antiinstantons, their interaction being not 
very strong (thus the name "liquid"). Instanton density is approximately 1 fm -4 . This model solves 
a number of problems in QCD, in particular it naturally leads to SCSB, explains the rf meson 
mass. However, there is a number of problems of fundamental importance: it remains unclear, 
how to make such instanton-antiinstanton ensemble stable (infrared inflation of instantons), and 
confinement cannot be explained in the instanton liquid model. 

However, apart from quasiclassical instantons there are other nonperturbative fields in the vac- 
uum, which allow to solve "infrared problem" of instantons. Nonperturbative vacuum of QCD 
can be parametrized by gauge invariant vacuum averages of gluon field strength tensor - vacuum 
correlators 3 . This allows to successfully describe many phenomena in QCD (see the review pa- 
pers 4,5). Bilocal correlator is sufficient for qualitative description of many phenomena (stochastic 
vacuum model). Moreover, there are indications that corrections due to higher correlators are small 
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and in some cases amount to several percent jlj. Nonzero string tension naturally appears in this 
approach and can be expressed through vacuum average of chromoelectric field [3] 



where Ei — Fa is chromoelectric field and F^ v = F^ v t a is the field strength tensor of gluon field. 
$(a;,0) is the parallel transporter along the straight line, connecting points x and 0. Integration 
in goes in space-time plane (i,4). 

The problem of infrared inflation of instantons is also naturally solved in the stochastic vac- 
uum model. Investigation of influence of nonperturbative quantum fluctuations on instantons was 
started in papers [SJ [7] . Later it was shown that standard perturbation theory changes in the 
nonperturbative stochastic vacuum, and contribution of large instantons to physical quantities 
becomes finite $ . Next, direct interaction of instanton with nonperturbative vacuum fields was 
shown [5] to lead to stabilization in instanton size, the average size of instantons being of the 
order of vacuum correlation length. Distribution of instantons over sizes was found to be in good 
agreement with lattice data, with a maximum at p c ~ 0.25 -j- 0.3 fm. 

An important property of vacuum is the scaling of string tensions between sources in different 
representations of color group, which is observed in the lattice calculations JH]- It was shown 
in that instanton contribution to the force between heavy quarks in different representations 
violates Casimir scaling, and this leads to a constraint on instanton density. 

In the present paper we study instantons' contribution to the bilocal correlator at finite tem- 
perature 1 . We demonstrate that correlation length of bilocal correlator, evaluated on calorons 
(calorons are instanton-like configurations at finite temperature) in the lowest order in caloron 
density has significant dependence on temperature. At the same time it is found in the lattice 
calculations that within errors correlation length for chromomagnetic fields does not depend on 
the temperature in the whole region between zero and critical temperature T c |14| . 

Temperature dependencies of the chromomagnetic bilocal correlator and of the spatial string 
tension, a s (T), were found analytically in the paper |15|. It was shown there that chromomagnetic 
condensate slowly grows with increasing temperature up to 2T C , while correlation length X m (T) 
does not depend on T in this temperature region. At high temperature, T > 2T C , condensate grows 
with T as {H 2 ) T oc g 8 (T)T 4 and magnetic correlation length decreases as A m (T) oc \/{g 2 (T)T). 
Obtained analytical expression for the cr s (T) agrees well with lattice results for all tempera- 
tures 16 . There exists a simple physical explanation of the independence of correlation length on 
the temperature at T < T c . Inverse correlation length M = 1/A m is a mass of lightest gluelump 
1^ . Gluelumps |17| are not physical objects, however they play an important role in the nonper- 
turbative QCD, because their masses determine field correlators, and in particular, string tension 
at T — is related to gluelumps. Their masses were calculated analytically in the framework of 
QCD sum rules [T^], QCD string model and were determined numerically in various lattice 
calculations [201 - Summing up all those results, correlation length equals X m « 0.15 fm, which 
corresponds to M « 1.5 GeV. Taking into account that phase transition temperature is much 
less than gluelump mass T c -C M, it is physically clear that temperature dependence of gluelump 
mass, and therefore of the correlation length, should be very weak in the low temperature region 
T < T c ~ 300 MeV, which is confirmed by analytical ^5] an d lattice calculations. It should be 
mentioned, that magnetic correlation length does not change it's behavior at the phase transition 

1 Bilocal correlator in the instanton gas at zero temperature was evaluated in I12II13I . 
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point. It is well known, that order parameter for the deconfincment phase transition is vacuum av- 
erage of Polyakov loop, and therefore it is Polyakov loops correlator (L(x)L^ (0)) that significantly 
changes it's behavior at temperature T = T c . 



2 Chromoelectric and chromomagnetic correlators 

At finite temperature euclidian 0(4) space-time symmetry is broken to space 0(3) symmetry, and 
bilocal correlator is described by independent chromomagnetic and chromoelectric correlation func- 
tions, which are expressed through four independent functions D E (x 2 ), Df(x 2 ), D B (x 2 ) Df(x 2 ): 

f f)D E \ BD E 

g 2 (Tr(E i (x)$(x,y)E j (y)&(x,y)))=S ij (L> e + Df + z 2 4 ^ J + 

g 2 (Tr(B i (x)$(x,y)B j (y)&(x,y))) = Sii [D B + Df + z 2 ^-f j - KZ^, 

where Bi = \sijkFjk is the chromomagnetic field, Ei = Fa - chromoelectric field and F^ = F^ v t a 
is the field strength tensor. 



$(x, y) = Pexp (^igz^ J ds A p 



(y + sz)) (3) 



is the parallel transporter along the straight line, connecting points x and y; z — x — y. These 
phase factors $ are introduced to ensure gauge invariance of correlators. Obviously Df = Df = 0, 
D E = D B for the selfdual field (see ^3]). Thus it is sufficient to consider chromomagnetic correlator 
only. 

The very fact that for instantons D E = D B implies that instantons cannot be the only nonper- 
turbative fields in the vacuum of gluodynamics, because it is known from lattice calculations that 
chromoelectric and chromomagnetic condensates have different dependences on temperature. In 
particular, at the deconfincment phase transition temperature chromoelectric condensate abruptly 
vanishes, while chromomagnetic condensate stays intact |14| . 

3 Calorons' contribution to the chromomagnetic correlator 

Instanton field with the center at the origin in the singular gauge is given by 



A ^(x) = --v%d v hxR^ st \x), 



(4) 



„2 J2 



n( in8t >(a;) = 1 + fL = 1 .i p 

where x = (r, r). This expression is easily generalized to the case of finite temperature, and caloron 
has a well known form [21| : 

A;(x) = -^fj;,d u \nU(x), 

u(x) i + V p2 - 1 ■ Kp2 sh(27rr//3) (5) 



n— — oo 



r 2 + (t - n/3) 2 fir ch(27rr//3) - cos(2ttt//3) 
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where (3 = 1/T is the inverse temperature. Caloron field is given by 



F a — r) A a — f) A a A- np abc A b A c 

X ' XX' 

9 B i = -5ai9i(x) +e ai j-^-g 2 {x) + -^-gs{x), 
K = B i, 

,. 1 1 an i d 2 u i fon\ 2 1 /an\ 2 

5i(a ; ) = -771Z: + T71^T-T72 h^r tt2 ( I ' ( 6 ) 



rU dr U dr 2 U 2 \8t J U 2 \dr ) ' 

2 an an 1 a 2 n 



92{x) = - = 



(7) 



II 2 dr <9t II drdr ' 

2 /<9n\ 2 ia 2 n 1 1 on 
93( - x > ~ W (jfr ) _ iiF + rn¥' 

Owing to the selfduality of caloron field F^ u — F^ v chromoelectric field is equal to chromomagnetic 
one. 

Next step is the evaluation of the parallel transporter on caloron: 

$(x,y) = Pexp (i9^ ds A^y + szfj = Pexp (-iz^J a ds d v \nU(y + sz)j = 

= Pexp (^-it a ds [fffjZidj Inn + ffazidi Inn + fftjZidj lnll]j , 
z = x — y. 

For comparison with lattice results we will consider the case Z4 = 0. Then 

y) = Pexp (^-it a ds [fjtjZidj In n + tf 4 z t d 4 Inn] ^ . (8) 

In general one has to evaluate path ordered integral. However there are special cases, when 
the exponent commutes at different values of s. Obviously, these cases are r = 0, t = (3/2 and 
t = (3, since cMnnfr, r = 0;(3/2;(3) = 0. Another example is the single instanton field, since 
9 4 lnn( inst ) = T(l/r)(dlnn( inst V<9r). 

To evaluate the contribution of caloron gas to the two-point correlator in the lowest order in 
density one has to average single-caloron contribution over it's position, or, which is the same, to 
average over y at a given z = x — y. 

Averaging over y comes to the integral J d A y = J d 3 y dy 4 . Averaging over three-dimensional 
vector y docs not make any problems, but the integral over y 4 is not tractable, since correlator can 
be numerically calculated only when r = 0,(3/2,(3. Nevertheless it will be seen from comparison 
with the single-instanton case that at low temperatures calculations with r = arc sufficient, 
and at high temperatures correlation lengths for t = and r = (3/2 cases coincide, which makes 
averaging trivial. 
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4 Numerical results 

Let us use the equation 

g 2 {Tr(B i (x)*(x,y)B j (y)*Hx,v)))i; = 8ij (* + di + 2°^) - z iZj ^, (9) 

at z 4 = to define the functions d and di, which depend on z 2 and j/4 = r: d = d(r, z 2 );e?i = 
di(r, z 2 ). Averaging over y\ yields 

AS 

/ drd(T,z 2 ) ^D B (z 2 ) 1 
Jo 

f dTd 1 (T,z 2 )=Df(z 2 ) = 0, (10) 
Jo 

4tt 2 

D>V = 0) = ^. 

Last equality follows from the relation g 2 j d 3 y dj/4 (F^ u ) 2 — 32ir 2 (the field is selfdual, and this 
integral is proportional to the topological charge ~ J d A yFF). 

Let us define correlation length A T as the coefficient in the exponent of the best fit to the 
function d(r, z 2 ): 

d(T,Z 2 ) ~e~ WK (11) 




0.2 0.4 0.6 0.8 1 1.2 1.4 



z (fm) 

Figure 1: d< inst ) (r, z 2 ) at r = (0,0.15,0.3,0.45,0.6) fm; p = 0.3 fm. Upper line corresponds to 
t = 0. 

As discussed above, we will evaluate d(r = 0, z 2 ) and d(r = /3/2,z 2 ) only. However, this 
information is sufficient to determine the correlation length of the function D B (z 2 ). Indeed, by 
averaging over r one can find that normalized to 1 functions d(r — 0, z 2 )/d(r = 0, z 2 = 0) and 
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Figure 2: Correlation length of functions d(r = 0, z 2 ) and d(r = (3/2, z 2 ): X T= o((3) (solid line) and 
A T= £/2(/3) (dashed line); p = 0.3 fin. 

D B (z 2 )/D B (z 2 = 0) practically coincide. The fact is that instanton is a well localized field config- 
uration, and it's field decreases very fast when moving away from it's center. This is illustrated 
in the Fig.^ It shows function d( inst ) (r, z 2 ), calculated for the single instanton configuration for 
several values of r. Profile of the function D B (z 2 ) is determined by the function d(r — 0, z 2 ). The 
same is true for caloron when (3 3> p, because in this case instantons forming caloron are well 
separated. At higher temperatures, when (3 < 2p the functions d(r = 0,z 2 ) and d(r — (3/2, z 2 ) 
have almost the same correlation length (see Fig.[5J), which makes averaging over r trivial (it only 
affects the amplitude of the functions, which can be found from the condition D B (z 2 = 0) = ^f-)- 

Taking into account all aforesaid, we can restrict our analysis by the function g?(t = 0, z 2 ) only. 
It's correlation length is close to the correlation length of function D B (z 2 ) for all values of (3. 

Chromomagnetic correlator for r = 0,(3/2,(3 takes the form 

g 2 (Tr (B i (xMx,y)B j (y)^(x,y))) i j = g 2 J d*y (B?{x)U a \x,y)B b 3 (y)) , 

U ab (x, y) = Tr (t a $(x, y)t b &(x, y)) = ^S ab cos(20) - \e abc n c sin(20) + n a n b sin 2 (0), (12) 

n a ^\n%z lV3 J Q ds (± dlnUi l + SZ) ) 

Correlator defined in (|12|l is presented in Fig. 01 for several values of (3. Corresponding correlation 
lengths are shown in Fig. 0] 

5 Discussion 

Bilocal correlator at finite temperature has been studied on lattice in the paper |14j. and it was 
found that bilocal chromomagnetic correlator practically does not depend on temperature in the 
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Figure 3: d{r = 0,z 2 )/d(r = 0, z 2 = 0) for /3 = (0.05, 0.1, 0.2, 0.3, 0.5, 0.86, oo) fm" 1 ; p = 0.3 fm. 
The largest correlation length (upper line) corresponds to zero temperature {(3 — oo), i.e. to 
instantons' contribution. 



whole region from zero to T c (this dependence is also very weak above critical temperature). At 
the same time, as it is seen from Fig. QJ the correlation length in the dilute instanton gas changes 
by approximately 20% when temperature changes from zero to 300 MeV. Thus, one comes to the 
conclusion that instanton gas cannot be a correct description of the real vacuum of gluodynamics. 
There are several possible explanations of these lattice data. 

1) Instanton density is much less than 1 fm -4 , and the main contribution to the correlator comes 
from other nonperturbative fields, which are responsible for the confinement. 

2) Instanton density is high enough to make next terms in density expansion change the result 
qualitatively leading to the independent of temperature correlation length. 

3) Instanton ensemble is rearranged at nonzero temperature, instantons and antiinstantons form 
molecules, and the contribution of these instanton-antiinstanton molecules to bilocal correlator is 
significantly different from that of instanton gas. This binding may be connected to the decon- 
fincmcnt phase transition. For example, this is the case for 3-d adjoint Higgs model, where it was 
found that in the deconfined phase instantons and antiinstantons are bound into molecules |22j . 
Formation of instanton-antiinstanton molecules in QCD vacuum and it's connection to the chiral 
symmetry restoration phase transition is discussed in review |23| . 

4) It should be noted, that in the last years a different scenario of the caloron vacuum in the 
confined phase (T < T c ) is discussed. There exist field configurations with nontrivial holonomy, 
the so called KvBLL solution [21]. This scenario suggests that calorons may be composites 
of dyons and antidyons. Lattice calculations, giving numerical evidence for the existence of these 
field configurations for the cases of SU(2) and SU(3) groups, are presented in papers |2lfl \Tf\. 
Clearly, in the case of this scenario correlation length in the nonperturbative vacuum has a more 
complicated temperature dependence, than in the case of caloron gas studied in this paper. 

In principle, one can consider one possibility more, when instanton size distribution depends on 
temperature in such a way, that the resulting correlation length does not depend on T. However, 
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Figure 4: Correlation length as a function of temperature A r= o(T); p — 0.3 fm. 



this is ruled out by existing lattice calculations, which demonstrate that within errors instanton 
size distribution stays the same in the temperature region from to T c |28| . 

In this paper we have found the contribution of instanton gas to the bilocal correlator at finite 
temperature. Comparison of the temperature dependence of correlation length with lattice results 
allows to make the conclusion that if instanton density is not too high and corrections to the dilute 
gas approximation are not significant, then instanton density should be much less than 1 fm~ 4 , 
the value standard for instanton liquid model. Vacuum of non-Abelian gauge theory at finite 
temperature is a very complicated system, and it is natural to expect that it can be a combination 
of the above listed possibilities. 

Authors are grateful to Yu.A. Simonov for useful comments. This work is supported by NSh- 
1774.2003.2 grant and by Federal Program of the Russian Ministry of Industry, Science and Tech- 
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